Active and Passive Fields in Turbulent Transport: the Role of Statistically Preserved 

Structures 
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We have recently proposed that the statistics of active fields (which affect the velocity field itself) 
in well-developed turbulence are also dominated by the Statistically Preserved Structures of auxiliary 
passive fields which are advected by the same velocity field. The Statistically Preserved Structures 
are eigenmodes of eigenvalue 1 of an appropriate propagator of the decaying (unforced) passive 
field, or equivalently, the zero modes of a related operator. In this paper we investigate further this 
surprising finding via two examples, one akin to turbulent convection in which the temperature is the 
active scalar, and the other akin to magneto-hydrodynamics in which the magnetic field is the active 
vector. In the first example, all the even correlation functions of the active and passive fields exhibit 
identical scaling behavior. The second example appears at first sight to be a counter-example: the 
statistical objects of the active and passive fields have entirely different scaling exponents. We 
demonstrate nevertheless that the Statistically Preserved Structures of the passive vector dominate 
again the statistics of the active field, except that due to a dynamical conservation law the amplitude 
of the leading zero mode cancels exactly. The active vector is then dominated by the sub-leading 
zero mode of the passive vector. Our work thus suggests that the statistical properties of active 
fields in turbulence can be understood with the same generality as those of passive fields. 

PACS numbers: 47.27.-i 



I. INTRODUCTION 

The aim of this paper is to address the statistical 
physics of so called "active" fields in developed fluid tur- 
bulence. These arc fields that differ from the fundamen- 
tal fluid velocity field u(r,t), but that interact with the 
velocity field in an essential way, for example effecting a 
significant change in the scaling exponents of the veloc- 
ity correlation functions from the classical Kolmogorov 
exponents. For the sake of concreteness we will focus 
on two generic examples with very different interactions 
between the active and the velocity fields. 

The first is thermal turbulent convection, in which 
the temperature field T(r, t) is driving the velocity field 
through buoyancy effects. In the Boussinesq approxi- 
mation the temperature equation reads like a standard 
forced scalar advection problem, 



changes the scaling exponents of u. When the conditions 
are right it may even change the scaling exponents from 
Kolmogorov to Bolgiano (up to anomalies) [y . 

The second example is that of magneto- hydrodynamics 
(MHD), in which the magnetic field b(r,t) is driving the 
velocity field u(r, t) according to 



-Vp + b-Vb + vV 2 u 



du 

— + u- Vu 

at 

db 

— +u-Vb = b-Vu + nV 2 b + f 
dt 



(3) 



dT(r,t) 
dt 



-u(r,t) ■ VT(r,t) = KS7 2 T(r,t) + f(r,t). (1) 



Here K is the thermal diffusivity and f(r,t) is a white 
random force of zero mean with compact support in k- 
space, acting on the largest scales of the order of the 
outer scale L only. The velocity field is affected by the 
temperature. For an incompressible fluid of unit density 
[PJ (dropping the dependence on (r,t) for brevity), 



These equations of motion conserve (in the inviscid, un- 
forced limit) three quadratic invariants, i. e. the energy, 
magnetic helicity and cross helicity || 

Our main interest is in the properties of the statisti- 
cal objects characterizing the active fields, including their 
anomalous scaling. Here "anomalous scaling" means that 
multi-point correlation functions are homogeneous func- 
tions of their arguments, with exponents that cannot be 
guessed from dimensional analysis. Thus for example the 
field <fi(r, t) (with <j> being T or b respectively) has simul- 
taneous multi-point correlation functions 



i) = (<t>(ri,t)<j>(r 2 ,t)- ■ -4>{r m ,t)) f 



du 
~dt 



u ■ Vit = — 'Vp 



vV 2 u 



agTz 



(2) 



Here p, v, a, g and z are the pressure, kinematic viscosity, 
volume expansion coefficient, acceleration due to gravity 
and a unit vector in the upward direction respectively. 
The appearance of T in the equation for u is crucial, and 



( 4 ) 

where pointed brackets with subscript / refer to averag- 
ing over the statistics of the advecting velocity field and 
of the forcing. The forcing is taken to be white random 
noise with zero mean. When the forcing is stationary in 
time this object is time independent. Anomalous scaling 
means that 



F (m \X 



ri, 







(5) 
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with £ m having a non-trivial dependence on m. In what 
follows we will assume that the advecting velocity field 
itself is fully turbulent, and that its correlation functions 
are also exhibiting scaling behavior like Eq. (||) . 

The main point of this paper is that the statistical the- 
ory of the active fields calls for consideration of auxiliary 
passive fields that satisfy the same equations of motion 
as the active fields, but do not affect the velocity field 
itself. In other words, For the two problems at hand we 
consider the following equations of motion : 



du 

~dt 

dT 

~dt 
dC_ 

~dt 



u ■ Vu 
u-VT 
u-VC 



on the one hand, and 

du _ 

— + u ■ Vu 

at 

— + u-Vb - 

ot 

m +u ' Vq - 



— Vp + vV 2 u + agT'z , 
«V 2 T + / , 

«V 2 C + / , (6) 

-Vp + b ■ Vb + i>V 2 u , 

b-Vu + nV 2 b+ f , 

q Vu + K \7 2 q + f , (7) 



on the other. 

Note that the velocity field that appears in the equa- 
tions for the passive fields is the same as the velocity field 
that results from solving the coupled equations of the as- 
sociated equations for the active fields. The forcing terms 
/ in Eq. (g) (resp. / in Eq. (Q)) have the same statistics 
as the forcing terms / in Eq. (|l|) (resp. / in Eq. (||)), 
but they must have different realizations. While it is not 
true of course that the statistics of the passive fields are 
independent of the statistics of the velocity fields, it is 
true that the statistics of the velocity fields are indepen- 
dent of the statistics of the passive forcing terms. This is, 
however, not the case with the active forcing terms since 
these forcing terms affect the active fields that affect in 
their turn the velocity fields. It is thus not at all evident 
at first sight that there should be any relation, apriori, 
between the statistics of the active fields and their pas- 
sive counterparts. On the other hand, if there were such 
a relationship, this would be very advantageous, since the 
statistics of the passive fields is understood as explained 
next. 

To understand the progress made in the context of pas- 
sive fields H| H , note that the passive fields satisfy a linear 
equation of motion that can be written as 



d<t>(r,t) 
dt 



= £4>( r> t) + f(r,t) 



(8) 



with the actual form of the operator C determined by the 
problem at hand. In recent work j^, Q it was clarified why 
and how passive fields exhibit anomalous scaling, when 
the velocity field is a generic turbulent field. The key is 
to consider a problem associated with Eq. (0) which is 



the decaying problem in which the forcing f(r,t) is put 
to zero. The problem becomes then a linear initial value 
problem, 



d(j)/dt = C<t> , 



with a formal solution 



tj)(r,t) = j dr'R(r,r',t)(f>{r',Q) 
with the operator 

R = T+ cxp[ / dsC(s)} , 
Jo 



(9) 



(10) 



(11) 



and T + being the time ordering operator. Define next 
the time dependent correlation functions of the decaying 
problem: 



G^(n,---,r m ,t) = {4>{r x ,t)- 



,*)). (12) 



Here pointed brackets without subscript / refer to the 
decaying object in which averaging is taken with respect 
to realizations of the velocity field only. As a result of Eq. 
( |io| ) the decaying correlation functions are developed by 

a propagator V^J , (with r = r 1 ,r 2 ,- ■ -,r m ) : 

G^ m \r u --;r m ,t) = JdpP ( ^(t) G(" l )( Pl ,-.-,p m ,0) . 

(13) 

In writing this equation we made explicit use of the fact 
that the initial distribution of the passive field 0(r,O) is 
statistically independent of the advecting velocity field. 
Thus the operator can be written explicitly 



r |p 



P$(t) = (R{r 1 ,p 1 ,t)R(r 2 ,p 2 ,t)---R(r m ,p m ,t)) 



(14) 



is that the operator V^" l J pos- 



The key finding | 

sesses left eigenfunctions of eigenvalue 1, i. e. there exist 
time-independent functions Z^ m \ri, r 2 , ■ ■ ■ , r m ) satisfy- 
ing 

Z^( ri , ■ ■ -,r m ) = J dpT^(t)Z< m H Pl , ■ ■ -,p m ) . 

(15) 

The functions Z' m ) are referred to as "Statistically Pre- 
served Structures" , being invariant to the dynamics, even 
though the operator is strongly time dependent and decay- 
ing. How to form, from these functions, infinitely many 
conserved variables in the decaying prob lem was shown 
in 0, and is discussed again in Sect. [II B| . The functions 
are homogeneous functions of their arguments, 
with anomalous scaling exponents £ m : 



(16) 



where ". . ." stand for subleading scaling terms. Since 
Eq. fll~5| ) contains Z^^r) on both sides, the scaling ex- 
ponent £ m cannot be determined from dimensional con- 
siderations, and it can be anomalous. More importantly, 



3 



it was shown that the correlation functions of the forced 
case, F( m )(r) Eq. (|4), have exactly the same scaling ex- 
'(21)0. 



ponents as 



i (r ) |7J] . In the scaling sense 



(17) 



This is how anomalous scaling in passive fields is under- 
stood. Lastly, we note that for the operator governing the 
time derivative of Eq. (|l|), Z^ m \r) is a zero mode. We 
will use the terms "Statistically Preserved Structures" 
and "zero modes" interchangeably. 

Of course, returning to the active fields, it makes no 
sense to consider the decaying problem; as the active field 
decays, the statistics of the velocity field changes, and 
there is very little to say. On the other hand, we propose 
that it is possible to learn a great deal from consider- 
ing the forced solutions, comparing the forced correlation 
functions of the active field with those of the passive field 
when advected by the same velocity field |q]. The rest 
of this paper is devoted to making this point clear and 
solid. 

In Sect. H we discuss the active problem (||) in compar- 
ison with the passive problem (|^). A preliminary report 
of the correspondence between these problems was pre- 
sented in ||] . Since we are interested in points of princi- 
ple rather than quantitative details, we opt to work with 
a shell model of the turbulent convection problem. We 
will argue (cf. Sect. IV) that there are excellent rea- 



sons to believe that the results found for the shell model 
translate verbatim to the partial differential equations. 
The main result of Sect. |l| is that the forced 2mth- 
order correlation functions of the active and passive fields 
are both dominated by the Statistically Preserved Struc- 
tures of the deca ying passive problem, i. e. the functions 
Z^ 2m \r) of Eq. (|15|). The anomalous scaling exponents 
are the same for the passive and active forced correlation 
functions, they are universal (independent of the forc- 
ing f(r,t)) and determined by the scaling exponents of 
Z^ 2m \r). We present a careful discussion of the role of 
the statistical correlations between the forcing and the 
velocity field that exist in the active case, but are absent 



in the passive case. In the present problem the net result 
of these correlations is just an amplitude factor relating 
the moments of the two fields. In Sect. Ill 
a shell model of magneto-hydrodynamics 



we turn to 
On the face 



of it, this is a counter-example to the previous case: the 
active and passive fields exhibit radically different scal- 
ing exponents. The main result of this section is that 
nevertheless the Statistically Preserved Structures of the 
passive problem are shown to dominate the statistics of 
the active problem, but the existence of a conservation 
law in the latter results in an exact cancellation of the 
amplitude of the leading zero mode. We identify analyti- 
cally the leading and subleading exponents of the passive 
problem, and then observe the cancellation of the leading 
contribution by the dynamics. The Summary section IV 
presents the general lesson for the statistical physics of 
the (nonlinear) active problem. We propose that the zero 
modes of the auxiliary passive fields will always have a 
dominant role in the statistics of active fields. The active 
fields will thus share the same scaling exponents as the 
passive fields unless there exist additional conservation 
laws for the active fields. In all cases the calculation of 
the active scaling exponents can be achieved in the con- 
text of the passive problem, which boils down to finding 
the zero modes of a linear operator. 



II. ACTIVE AND PASSIVE SCALARS IN A 
MODEL OF TURBULENT CONVECTION 



A. Model and numerical results 

In this section we examine in detail a shell model of 
active and passive scalars for which the statistical object 
can be computed to high accuracy. We consider a model 
that reproduces the conservation laws and the form of 
coupling between the active field and the velocity field 
in Eqs. (Q) and (0). Our model is a variant of the shell 



model studied in ref. |10 



= ak n {u 2 n ^ x - \u n u n+ i) + 6fc„(u„u„_i - Xu 2 l+1 ) - vk 2 n u n + T n , (18) 
= M n (u n -xT n _x - \u n T n+1 ) + bk n (u n T n ^x - Xu n+ xT n+1 ) — Kk n T n + f 6„ t o , (19) 

dC 

= afc„(u„_ 1 C„_i - Xu n C n+1 ) + 6fc„(u„C„_i - Xu n+1 C n+1 ) - nk 2 n C n -I- f S n ,o ■ (20) 

I 



In this model all the field variables are real and n stands 
for the index of a shell of wavevector k n = ko\ n , with 
n = 0, 1, ■ • • , N — 1. We take A = 2, and the parameters 
used in the simulation are a = 0.01, 5 = 6 = 6=1, 
feo = 1, k = v = 5 x 10~ 4 . The number of shells is 



N = 30, and the forcing is white noise of zero mean on 
the first shell. 

Without the coupling to T n , the velocity equation has 

— 1/3 

an inviscid unstable Kolmogorov fixed point, u n ~ k n 



4 




n 



FIG. 1: The scaling exponents £JJ of the velocity field (circles) 
and of the active scalar field (squares) for even n's. The 
solid lines are respectively 3n/5 and n/5 for the velocity and 
the active scalar fields. Shown in the inset is the PDF of 
z = u n /{ul} 1/2 at shell n = 14. 

This is changed by the coupling pc] , and the system of 
equations for T n and u n exhibits an inviscid unstable 
Bolgiano fixed point, u n ~ fc„ 3 ^ 5 , T n ~ k n X ^ . The 
chaotic dynamics renders the statistics of the velocity 
field strongly non-Gaussian (cf. inset in Fig. [j]). The 
exponents for the active scalar are markedly anoma- 
lous, whereas, for the velocity, they appear closer to nor- 
mal (see Fig. [I]). The equation of motion for the passive 
field C is identical to the equation of motion of T, but 
it does not affect the velocity field u. This equation has 
a C — > — C symmetry, whereas the coupled system of T 
and u lacks this symmetry. This difference is reflected in 
the statistics of the two fields. 

To demonstrate this difference between the active and 
passive fields, we show in Fig. || the probability distri- 
bution functions (PDF) of x = 4> n / i&n) 1 ^ 2 where <pn is 
T n or C„, for n = 14. One clearly sees the symmetry 
of the PDF of the passive scalar, in contradistinction to 
the asymmetry of the PDF of the active scalar. This is 
typical to all n in the inertial range. This is a demon- 
stration of the discussion after Eq. ( p8| ) . For the passive 
scalar the odd moments vanish, whereas for the active 
scalar they all exist. The situation is altogether different 
for the statistics of even moments. To demonstrate the 
difference we plot in Fig. | the (typical) PDF of f% and 
C 2 n for n = 9 and 14, where T n = T n - (T„). In plotting 
we realize that the passive scalar is defined up to a con- 
stant, so for the passive scalar the PDF is plotted for the 
rescaled variable /3C^, where 

^ = (T„ 2 )//(Cn)/ ~ 0.6327 . (21) 








-15 -10 -5 5 10 15 

x 

FIG. 2: The PDF's of the active (solid) and passive (dashed) 
scalars at shell n = 14. Note that the PDF of the active scalar 
is asymmetric. 
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-1 1 1 1 1 1 ■ 1 1 1 1 1 

1 2 3 4 5 

10 4 j 

FIG. 3: The PDF's of y where y = f £ (solid) or (iCl 
(dashed) at shells n — 9 and 14. 

Note that there is only one numerical freedom /3, constant 
for all n in the inertial range. An understanding of this 
numerical constant based on dynamical considerations is 
given in the next subsection. 

We find very close agreement of all the PDF's in the 
inertial range. The identity of the PDF's of T% and C% 
translates automatically to the identity of the even-order 
structure functions F^ 2m \k n ) = (^ m )/, where 4> n = T n 
or C n (up to a constant f3 m ). This is demonstrated in 
Fig. |J. We see that the 2nd, 4th and 6th-order structure 
functions are barely distinguishable, with the same scal- 
ing exponents in the inertial range. Finally, we demon- 
strate that the identity of the statistics of the squares of 
the passive and active scalars transcends structure func- 
tions. Consider for example the multi-point correlation 
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FIG. 4: The even-order structure functions (T^ m ) / (circles) 
and (/3 m C^ m )/ (squares), with m = 1, 2 and 3, from top to 
bottom. 




FIG. 5: Upper: log 10 (T^ +6 )/ (circles) and 
logio^CnCn+s)/ (squares). Lower: log 10 {f^f^ +5 fi n ) f 
(circles) and log 10 (/3 3 C^ C^ +5 C| n )/ (squares). 



functions (T^ 



s ) and (f 2 f 2 +5 f 2 n ). 



^ \- n -n+ 5 -2«/- In Fi g^ 1 these 
correlation functions are compared to their passive coun- 
terparts. The conclusion is that again the multi-point 
correlation functions are indistinguishable once the pas- 
sive ones are rescaled by (3 q where q is the overall order 
of the correlation function. 



B. Analysis of the results 

To understand the results we start with the passive 
field, demonstrating that its forced structure functions 
are actually Statistically Preserved Structures. Consider 
then the decaying passive problem, i. e. Eq. (^0|) without 
the forcing term. The initial value problem for the time 



dependent structure functions G^ 2 " 1 ^ (k n ,t) = (C 2m (t)) is 
the shell analog of Eq. ( |l3| ) , 



G {2m \k ni t) = V^{t)G {2m \k n ,,t = 0) 



(22) 



which defines the 2mth-order propagator V^™) (*) • Here 
and below, repeated indices are being summed over. In 
Fig. H we show typical decay plots for the quantity 
K (2m) = Y,n G(2m H k n,t) for to = 1,2,3, starting from 
the initial conditions G' 2m ^(fc„, t = 0) = <5n,i6- 

Statistically Preserved Structures in this case represent 
left-eigenfunctions Z^ 2m \k n ) of eigenvalue 1 satisfying 



z( 2m \k n/ ) = Z^(k n )V^(t) 



(23) 



The statement that we want to demonstrate is that the 
forced structure functions F^ 2 " 1 ^ of the passive scalar 
scale like these eigen-modes of the decaying problem: 

(k n ) = (C 2m ) f ~ (k n ) . (24) 

To demonstrate this we use the method of || and define 
the quantities I^ 2m \ 



j(2m) = J2G (2m) {k n ,t)F < - 2r ^(k n ) 



(25) 



Using Eqs. (g2|) and (|23) we see that if Eq. ( p4|) is obeyed, 
than the quantities i^ 2 ™ 1 ' are time independent. Indeed, 
in Fig. H we demonstrate the stationarity of these ob- 
jects, thus supporting Eq. (]24|). The analytic explana- 
tion as to why the forced solutions agree with the Sta- 
tistically Preserved Structures of the decaying problem 
was provided in JtT]. Before turning to the active field, it 
is worthwhile to observe how any initial condition of the 
decaying passive field lands on the scaling solution that 
is represented by the Statistically Preserved Structure. 
Consider the initial value experiment that is reported in 
Fig. ^. Here we start, as an example, from the initial 
value C n (t = 0) oc kf/ 3 . In this initial condition the 
order of the amplitudes is inverted with respect to the 
spectrum of the passive scalar. We plot, as a function 
of time, the trajectories of C n (t) as computed just from 
this initial condition, averaged over 650 realizations. We 
see that the trajectories land on a decaying scaling solu- 
tion in which the order of the amplitudes and the ratios 
between them are identical to the spectrum of the zero 
mode of the passive field; the decay that we see, at a rate 
proportional to 1/i 2 , is entirely due to dissipative effects, 
as explained in some detail in |t]]. 

Finally, we need to understand how the forced active 
scalar T falls on the Statistically Preserved Structure of 
the decaying passive problem, and what is the origin of 
the factor (3 in (|2l]). To this aim we note that both 
equations for passive and active fields can be written as 



at 



c 



+ M 



II. 



(26) 
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0.5 1 1.5 2 2.5 3 3.5 4 4.5 
t 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 
t 




0.5 1 1.5 2 2.5 3 3.5 4 4.5 
t 



FIG. 6: The decaying objects K {2m) (solid lines) and the con- 
served objects j' 2m ) (dashed lines) as a function of time, for 
m — 1,2 and 3. Time is measured here in units of the largest 
scale eddy turn over time r = [fcox/S^fco)] -1 ~ 22. In panels 
b and c we include in dotted lines the quantity /' 2m ' in which 
we replaced F { - 2m) by its dimensional prediction [F (2) ] m . We 
see that using the dimensional exponent does not make I^ 2m > 
time invariant. 



This equation has a formal solution in the form 

Mt) = Rn,n<{t\$)<t>n>{t = 0) + / dTR n . (t\T) f {T) , 

JO 

(27) 

where R n ,n> (t\r) is the shell analog of the operator (pd]), 

Rn, n -(t\T) = (t+ f exp[£(s)]ds] . (28) 

The first difference between the active and passive fields 
is encountered when we take the average of this equa- 
tion. For the passive case, the average can be taken by 
decorrelating /o and R n , n ' ■ Since the mean of the force 
/o vanishes, we get 

(Cn)f = (29) 

Such a decorrelation is, however, not allowed in the active 
case since the forcing /o is correlated with T, which is 
itself correlated with u and thus with R n , n '- Hence 

{T n )f= (R n .n>(t\0)Tn>(t = 0))/+ / d T {R nfi {t\ T ) fo{ T )) f 

Jo 

(30) 

and T has a non-zero mean. Similarly, the passive scalar 
has zero odd moments and its PDF is symmetric. On 
the other hand, the active scalar has nonvanishing odd 
moments and its PDF is asymmetric. 

In spite of this great difference between the active and 
passive scalars, there is a close affinity between the ac- 
tive field and the Statistically Preserved Structures of the 
passive field. To see this, we note that the first term on 
the RHS of Eq. (^FJ) represents a decaying field. We ex- 
pect that if the initial condition T n (t = 0) has any com- 
ponent on the Statistically Preserved Structure of the 
passive field, it will quickly relax everything else and will 
land exactly on that solution. In this respect it is just 
the same as the initial value experiment reported in Fig. 

. In terms of the relative amplitudes of the different n 
shells there is nothing in the fate of the initial value term 
to distinguish the active and the passive fields. The sec- 
ond term on the RHS of Eq. (^FJ) is more subtle. First, 
we note that for every value of r we again face a de- 
caying experiment that takes place between the times r 
and t. In the language of the passive field, the integrand 
can be read from a decaying field with initial condition 
C n (t = t) = fo{T)S n .o. Indeed, in our simulations below 
this is precisely how we evaluate integrals of this type. 
We break the interval [0,t] into N sub- intervals {Ti}fL 1 , 
Ti = (i/N)t, and start a decaying experiment with initial 
conditions f{Ti)8 n fi. Measuring C n (t) and summing up 
all the contributions yields an approximation to the inte- 
gral. Every term in the integrand is expected to land, for 
most of the time t — r, on the scaling solution of the pas- 
sive field, in much of the same way that the initial value 
term does. Thus both terms in Eq. ( |30| ) are expected to 
scale like the passive field, which nevertheless itself has 
zero amplitude due to the symmetry. 
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50 100 150 200 250 300 350 400 450 
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FIG. 7: An example of the fate of an initial value term as a 
function of time, in units of tq. The initial amplitudes are 
inverted in order compared to the 2nd order zero modes of 
the passive field. Shown are shells n — 5, 7, . . . , 19. 




5 10 15 20 

t-s 



FIG. 8: Comparison of the integrands in Eqs (|32j) (the passive 
case) and ( |33[ ) (the active case) for n = 10. The plots are 
indistinguishable . 



The correlation effects that play a role for the active 
scalar will be responsible for the factor (3 that we discov- 
ered in the numerics. To see this we need to consider the 
2nd order structure functions: 



F^(k n ) = <[i?„,n'(i|0)<M* = 0)] 2 )/ 

+2(R n , n ,(t\0)q> n ,(t = Q) [ dTRn,o(t\T)f (T)) f 



(31) 



+ / dT'dr''(R n , (t\T')R n , (t\T'')f (T')fo(T'))f . 

Jo 

For sufficiently long time the first two terms, denoted be- 
low as X\ and Z2 respectively, do not contribute to the 
structure functions, and any difference between the active 
and passive fields must be ascribed to the last term. The 
last term, denoted here as Z3 , has a "diagonal" contribu- 
tion, which is obtained for r' = r" and an "off-diagonal" 
contribution, which is the rest of the integral for which 
t' ^ t" . For the passive field, R n , n ' and /o decouple 
and only the diagonal part exists. For the active field 
there is no decoupling. Denoting this term l3 )( z, it reads 
respectively 

I3, d = / ds{R nfi (t\s)R nfi (t\s)) f Jo (passive) , (32) 



Z3, d = / ds{R nfi (t\s)R n . (t\s)fo(s)f (s))f (active)(33) 
Jo 

In Fig. H we compare the integrands of these two ex- 
pressions, measured directly in our simulation as ex- 
plained above. We can see that there is not much 
difference between them; the diagonal term cannot be 
blamed for the factor (3. On the other hand, in Fig. ^ we 
show the full integral Z3 and compare it with its diago- 
nal term. We see that in the passive case the diagonal 
part is everything, whereas in the active case there is a 





FIG. 9: The integral X3 in comparison to the diagonal part 
Xz,d for n = 10. Upper panel: the passive field. The integral 
agrees with the diagonal part at all times. Lower panel: the 
active field. The deviations are due to the non- vanishing con- 
tribution of the off-diagonal integral, which is also displayed 
in the next figure. The dashed line in both panels represents 
the stationary value of the corresponding second order struc- 
ture function. 



difference. Lastly, we show that this difference is pre- 
cisely the source of the factor (3. In Fig. 10 we show 
([^(t) - Rn,n'{t\0)<Pn>{t = 0)] 2 ) - I 3 ,d for the passive 
and the active fields. The former fluctuates all the time 



around zero. The latter is positive initially, and then be- 
comes negative. For later times it saturates at a negative 
value that is precisely responsible for the factor (3. 

In summary, we find that the even correlation functions 
of the active and passive scalars share the same scaling 
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x 10 




FIG. 10: The off-diagonal integral for n = 10, computed as 
([(/>n(t) — Rn in i(t\0)<f) n >(t = 0)] 2 ) — 1z,d for the passive and 
active fields respectively. For the passive field (dotted line) it 
fluctuates around zero, while for the active field it begins pos- 
itive, and then turns negative. For longer times it saturates 
at a constant negative value, giving rise to the factor j3. 



exponents simply because the zero modes of the decaying 
passive problem dominate the statistics of both fields. It 
is thus possible to understand the anomalous statistics 



of the active field in the same way as that of the passive 
field. We believe that this is a significant result that 
should be put to further experimental and numerical tests 
in the PDE version of the problem. We will return to this 
point in the discussion. 



III. ACTIVE AND PASSIVE FIELDS IN A 
MODEL OF MAGNETOHYDRO DYNAMICS 



A. Model and Numerical Results 

In this section we examine a shell model that repro- 
duces the type of coupling and the conservation laws of 
Eqs. (||). We need to be careful about the dynamo ef- 
fect which we want to avoid in order to have station- 
ary statistics. We thus construct the model to mimic 
2-dimensional MHD, in which there is an inverse cas- 
cade of energy. Accordingly we need to have large scale 
damping terms in the velocity equation, and a force at 
intermediate scales. In all our simulations below we force 
both the velocity and the active fields on shells 10 and 11 
(denoted by u/), using white noise of zero mean. We run 
the model with 35 shells. The equations are an adapta- 
tion of the MHD shell model of Ol to the Sabra shell 



model 1 12 . All field variables are complex numbers: 



du 

—^j- = ik n [a\(u* n+1 u n+2 - b* n+l b n+2 ) + b(u*_ x u n+ i - b*_ x 6 n+1 ) - cA -1 {u n - 2 u n -i - b n _ 2 &n-i) 

+fn&n,n f + vk 2 n U n + vk~ A Un , 
— ^ = ik n [a\(u* n+1 b n+ 2 - b* n+1 U n+2 ) + b(u* n _ 1 b n+ l - b*_ ]«„+].) ~ cA~* (« n _ 2 &n-l - &n-2«n-l) 

+ fn6 n,n,f ~l~ ftk n b n + Hk n b n . 



(34) 
(35) 



r 



The coefficients a, b, c, a, b and c can be parametrized To mimic the magnetic helicity, we can write down a 
as follows : 



o=l, b =-S, c=-(l-S), 

a = l — 5 — S m , b = 5 m , c = 1 - S m . 



(36) 



This choice ensures the conservation of the total energy 
and "cross helicity" in the inviscid limit v = v = k = 
k = f = f = 0, 



n 

K = ]>>(<&„) . 



generalized quantity 



(39) 



with a > a fixed parameter. We demand conservation 
of this generalized "magnetic helicity" , together with ab- 
sence of dynamo effect. This implies 



5 > 1 



5 = l + \~ a , 
S m - -1/(A« - 1) 



(40) 



(37) On the other hand, when 5 < 1 one can have dynamo, 
and therefore no stationary statistics. 

(38) I n addition to the conservation laws the equations of 
motion remain invariant to the phase transformations 
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u n —> % exp(z0„) and b n — > b n exp(iip n ). The conditions ip n + ip n +i - <fi n +2 = . (44) 
are 

<t> n + (j) n+ i - <t> n +2 = , (41) This implies ip n = cf> n Vn. 

4>n + 4>n+i — V'n+2 = , (42) -j>he p ass i ve field is denoted by q n , whose evolution is 

Vfo + ^n+i — V'n+2 = , (43) given by an equation similar to Eq. (pq), i. c. 



dq„ 



= ik n [a\{u* n+1 q n +2 - <7 n+ iW n+2 ) + b{Un_ x q n +\ - qn-iVn+l) - cA 1 (u n - 2 q n -i ~ q n -2U n -i 



dt 

+ fnS n 



(45) 



40 



30 




FIG. 11: Structure functions of order 2 (x), 3 (+) and 4 (o) 
for the passive field (dotted line), active field (dashed line) 
and velocity field (solid line). The two vertical lines denote 
the forcing shells. Note that the scaling exponents of the 
active field and the velocity field coincide. The parameters 
are N = 35, a = 2, k = 0.0625, v = 10~ 12 , v = 10~ 3 . The 
forcing is white noise on shells 10,11. 



(|u ra | p } ~ k n , and the spectrum decays like a power law 
in the "inertial range" which is between the forcing and 
the dissipative scales. We estimate from the numerics 
f]2 = C2 ~ -67, 773 = £3 w 1.0 and 774 = (4 1.33, in close 
correspondence with the Kolmogorov dimensional predic- 
tions. The passive field has exponents, defined similarly 
by (|<?n | p ) ~ k n l3p , that are with a different sign! Its spec- 
trum is an increasing function of k n in the inertial range. 
We measure ~ —1-33, (3% ~ —2 and fa ~ —2.67. If we 
assume that the passive field lands on the Statistically 
Preserved Structures of the passive decaying problem, 
then it appears that the active field does not do so. 

In the rest of this section we will show that this is 
actually not a counterexample to the proposition that the 
active field lands on Statistically Preserved Structures 
of the decaying passive field. It does. What happens 
here is that, due to the conservation law Eq. (|3^), the 
amplitude of the leading Statistically Preserved Structure 
with the negative scaling exponent is exactly zero. The 
active field then lands on a sub-leading zero mode, which 
has standard, positive scaling exponent (the positive sign 
refers to r-space representation, as in Eq. 

B. Analysis of the results 



The fields q n and b n are advected by the same velocity 
field, however b n is active, while q n is passive. The in- 
viscid passive equation has only one conserved variable, 
i. e. Eq. ( ^9|) with q n replacing b n . It also satisfies the 
same phase relations as the active field. We want to know 
whether the scaling properties of b n are determined once 
again by the Statistically Preserved Structures of the de- 
caying problem of the passive field q n . 



In Fig. [lljwe show the spectra of the passive and active 
fields respectively, obtained from a direct numerical sim- 
ulation with the parameters as detailed in the figure leg- 
end. This appears like a striking counter-example to the 
results of the previous section: the two fields have totally 
different scaling behaviors. The active field has "stan- 
dard" scaling exponents T] p , defined by (|& n | p ) ~ k n r ' p , 
that coincide with those of the velocity field, defined by 



To gain insight into this interesting situation we note 
that the analog of Eq. ( |28| ) describes the dynamics of our 
active field b n : 



b n {t) = Rn tn ,(t\0)b n ,(t = 0) + / drR n , nf (t\r)f nf (r) , 

Jo 

(46) 

with an obvious re-definition of the present operator 
Rn,n> ■ It is very revealing to examine the time de- 
pendence of the two terms on the RHS of this equa- 
tion. We measure time in units of the eddy turn over 
time of the f orcing shell 10. This is defined as t\q = 
[Wd^iol 2 )]- 1 » 3.35. We will examine a forced sys- 
tem which began running at t = — oo, denoting a generic 
time as t = 0. In Fig. |lj panel a we show the time- 
dependence of the first term for 6 values of n in the in- 
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ertial interval. We see that the initial conditions rep- 
resent, as expected, a 'standard' spectrum in which the 
amplitude b n decreases as a function of n. As time pro- 
ceeds, the decaying term cannot recognize its being 'ac- 
tive' from being 'passive', and it switches rapidly to the 
Statistically Preserved Structure of the decaying passive 
field, characterized by a negative exponent. If it were 
not for the second term on the RHS of Eq. (ft6[), then 
b n would have landed on the same solution as q n . What 
about the second term? In panel b of Fig. [12] we show 
the n dependence of the term at time t — 3 x 10~ 4 . 
We see that also this term agrees, in its n dependence, 
with the negative exponent of the passive field. Yet, the 
LHS b n (i) fluctuates around decreasing amplitudes as n 
increases, meaning that the leading (negative) exponent 
exactly cancels between the two terms on the RHS of Eq. 
(|h|) . We demonstrate this cancellation in Fig. [l3| There 
we plot the real parts of the initial value term and the 
integral term at time t = 0.3. We see that the two terms 
cancel each other. The imaginary parts exhibit the same 
behavior. 

Next we need to understand this cancellation from the 
analysis of the equations of motion. With this analysis we 
will also show that the solution on which b n (t) is landing 
is also a Statistically Preserved Structure of the decaying 
passive field, albeit with a sub-leading scaling exponent. 

C. Statistically Preserved Structures of the passive 

field 

In the next subsection we will show that the velocity 
field attains a scaling solution with £3 = 1: 

S 3 (k„) = 3 < u„_iu„u,* +1 >y~ k~ x . (47) 

where denotes the imaginary part. In this section we 
will assume this, and examine what are the scaling solu- 
tions that agree with the existence of a 2nd order Statis- 
tically Preserved Structure for the passive field. We are 
not going to compute the anomalous scaling exponent 
exactly, but rather obtain their dimensional estimates. 
Since we are after a radically different apparent behav- 
ior, small numerical corrections are not our main concern. 
With this caveat in mind, we can calculate the exponent 
/?3 characterizing the third order structure function 

P 3 (k n ) = $t(qn-iq n q* + i) ~ • (48) 

The condition for the existence of the 2nd order Statisti- 
cally Preserved Structure is 

!<|g„| 2 )=0, Vn, (49) 

in the inviscid limit. Using Eq. J45| ) this condition gen- 
erates a number of third order quantities that need to be 
analyzed first. Denote therefore 

Q3,l(K) = $ < Wn-l9n9n+l > , 
Q 3 A k n) = 3 < q n -iu n q* n+l > , 
Q 3 , 3 ( k n) = 3 < q n -iq n u* + i > ■ (50) 




0.05 0.1 0.15 0.2 0.25 0.3 

t 




10 15 20 25 30 

n 



FIG. 12: Panel a: The fate of the modulus of the initial value 
term, averaged over 2, 000 realizations. Shown are shells 15, 
17, 19, 21, 23 and 25 (top to bottom from the left-most side). 
At time t = their relative amplitudes agree with the scaling 
exponent of the active field. As time progresses the decaying 
field switches to the relative amplitudes which agree with the 
scaling exponent of the passive field. Panel b: The modulus 
of four realizations of the integral as a function of n, for time 
t — 3 x 10 -4 fin unit of rio). Note that both terms of the 
RHS of Eq. (Uq) exhibit the same leading scaling behavior. 
This is canceled exactly as is demonstrated in Fig. ^ 

In order to construct scaling solutions for these objects, 
Dimensional consideration imply that the fields involved 
in the averages above have scalings u n oc A~™/ 3 and q n cx 
A - ' 33 ™' 3 . We infer the expressions 

QsAK) = \q \ 2 \u \k-W°+W\-^-iy 3 , 

QsAkn) = M 2 KI^- (2/33+1)/3 , 

QsAkn) = M>o|^ (2/33+lW3 A (/33 " 1)/3 • (51) 
We can thus rewrite 

QaAkn) = \- {p3 - 1)/3 Q 3 {kn) , (52) 
QsAkn) = Q 3 (k n ) , (53) 
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!3,3 



(k n ) = \^-^ 3 Q 3 (k n ) 



where Qs{k n ) scales like 



Q 3 (k n ) = A-( 2 ^+D«/3 . 



(54) 



(55) 



FIG. 13: the real part of the initial value term (dashed line) 
and integral term (solid line) as a function of k n . This is a 
demonstration of the cancellation of the leading order term 
in favor of the subleading one. The imaginary parts behave 
in the same way. 



J 



having these definitions in mind we derive, by demanding 
Eq. ( f49| ) and substituting the scaling form of Q3 Eq. (p35|), 
the equation 



5A -(2/3 3 -2)/3 (1 _ A (/3 3 -D/3 ) + g (A -Cfc-l)/3 _ A (fc-l)/3) + g A (2/S s -a)/3( A -Cft,-l)/3 _ 1) = Q . (56) 

I 



This is a fourth order polynomial in A^ 3-1 )' 3 . The four exponents fa of the order of (neglecting anomalies) fa — 
roots are 2/3, —4/3, —10/3. To know which of these is physical, we 

need to consider the fluxes supported by these solutions. 

The only flux that is relevant for the passive field is 
the magnetic helicity. For the case considered here with 
5 > 1 it can be conveniently computed at the shell M by 
Here three of the roots correspond to a priori physical evaluating 
solutions: 



A (ft>-i)/3 




(57) 



fa = { 1 - 3a/2 , 
1 — 3a . 



(58) 



<S> H 



2 dt \ fc« 



(59) 



In our simulations with a = 2 these results are fa = 1, Using the equations of motion to evaluate this object we 
-2 and -5 respectively. This is in agreement with spectral find 

I 



k M+l(^ < lMUM+\q*M + 2 > -3 < qMqM+lU* M+2 > ) 



+ k M "( 3 < qM-lU M q*M+i > +S < u M _iq M q* M+1 > ) 

I 



(60) 



We can evaluate now the magnetic helicity flux for the 
three scaling solutions (p8|). We find 



(61) 




We conclude that the third solution is unphysical, since it 
supports a flux that diverges with M. The first two solu- 
tions are allowed. With fa — —2 we get a constant flux; 
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this is the leading scaling solution, and is indeed realized 
in the simulations. The solution P 3 — 1 is subleading, it 
is associated with a decaying flux, and is asymptotically 
allowed. It is not observed in the passive field simulations 
simply because it is subleading. 

Our main point will be that the active field will in 
its turn land on the subleading Statistically Preserved 
Structure because the additional conservation laws ex- 
clude the leading one. We demonstrate this phenomenon 
in the next subsection. 



D. Why does the active field fall on a subleading 
zero mode? 



If we accept the general philosophy that active fields 
exhibit scaling behaviors that are determined by the zero 
modes of the auxiliary passive fields, then we should ex- 
plain here why in the present case the active field avoids 
the leading zero mode, and appears to land on the sub- 



J 



leading one. The answer is hidden of course in the con- 
servation laws, as we expose now. 

We first repeat the analysis performed in the previous 
section to find the consequence of the equation 



or, equivalently, 



dt 



< \u n \ 2 >=0 



(62) 



(63) 



Using now the definitions 

B 3 (k n ) = Sy(6„_i6„6* +1 ) ~ k~ Vs 

Q 3 (k n ) ee fc-(^3+Cs)/3 , ( 64) 

we obtain an equation that is analogous to Eq. (|56|), 



|^l-(2'73+C3)/3( 1 _ ^(V3-( 3 )/3\ + 5( / \-(»)3-C 3 )/3 ^0? 3 -C 3 )/3-) + g^-l+(2»73+f3)/3M-073-C3)/3 _ J) _ Q (65) 

I 



This is a fourth degree polynomial for A 173 / 3 if £3 is known. 
Obviously, if we simply substituted here (3 = 1 we would 
get the same predictions for 773 as obtained for (3 3 in Eq. 
(p8|). However, we have in this case an important addi- 
tional constraint that is absent in the case of the passive 
field, which can be inferred from the additional conser- 
vation equation 

d „ , 

(66) 



-3? < u* n b n >-- 
dt 







Repeating the analysis as above, and introducing a new 
object A 3 (k n ) = fc^ 7?3+2 ^ 3 ' yields the two equations 

a\B 3 {k n+1 ) + bB 3 {k n ) + cA- 1 B 3 (fc n -i) - , (67) 
a\ l -^-i3)/*Az{k n+l ) +bA 3 (k n ) 

+ cA -l + (r, 3 -C3)/3 A3(fcji _ i ) = ; (6g ) 

Solving this system together with Eq. (|6^) yields the 
scaling exponents 



C3, V3 



1 



1 + log A (a/c) 



(69) 



It is easy to check that, among the four possible combi- 
nations of this equation, the only solution allowed by Eq. 



(3 



V3 



1 . 



(70) 



We thus conclude that as far as the active field is consid- 
ered, the additional conservation law rules out the leading 
zero mode of the passive problem, leaving us only with 
the subleading mode which is observed in the simulations. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we considered the correspondence be- 
tween the statistics of active fields that are advected by 
a turbulent velocity field, and the statistics of an auxil- 
iary passive field that is advected by the same velocity, 
but does not affect it. The two examples were akin to 
turbulent convection and to magneto-hydrodynamics re- 
spectively. In the first example the conserved variables 
for the equations of the passive and active fields are the 
same. For the second example the active problem ex- 
hibits additional conservation laws. This was shown to 
be very significant in determining the respective statisti- 
cal physics of the two problems. 

The two examples appear very different in superficial 
examination. In the first example the even-order statis- 
tics of the passive and active fields turned out to be the 
same up to a single multiplicative factor (3, common to 
all orders. The forced structure functions of the active 
field scale with exactly the same exponents as the pas- 
sive field, which in turn are dominated by the leading 
zero modes of the decaying problem. We analyzed in de- 
tail the source of the multiplicative factor (3 and showed 
that it stems from the additional correlation effects be- 
tween the forcing and the velocity field that are absent in 
the passive case. Nevertheless these correlation effects do 
not cause a change in the scaling exponents. The gen- 
eral lesson that we would propose on the basis of this 
example is that whenever there exist a problem in which 
the equation of motion of the active field does not sat- 
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isfy additional conservation laws compared to the passive 
case, the former field will exhibit structure functions that 
are dominated by the leading zero modes of the latter. 
This point is also pertinent to the second example. Here 
the active equations possess additional conservation laws, 
and indeed the active and passive fields exhibit differ- 
ent scaling exponents. Nevertheless we argued that the 
structure functions of the active field are still dominated 
by the zero modes of the passive problem, but not the 
leading ones. The additional conservation laws results 
in exact cancellations in the contributions of the leading 
zero modes, and the active problem lands on the next 
allowed sub-leading zero mode of the passive problem. 

As a generalization, consider then a sufficiently turbu- 
lent velocity field which advects an active field, scalar or 
vector, which in its turn is forced by a force having a 
compact support in k space. An auxiliary passive field 
which is advected by the same velocity field can be em- 
ployed to find the zero modes of the operator involved 
in the passive decay problem. On the basis of the intu- 
ition gained with the examples presented above, we offer 
the following tentative conjecture: the forced structure 
function of the active field will exhibit scaling exponents 
that are the scaling exponents of the aforementioned zero 
modes. Whenever the conservation laws of the active and 
passive problems coincide, these will be the exponents of 
the leading zero modes. When the active problem has ad- 
ditional conservation laws, these will be the next-leading 
zero modes, as allowed by the conservation laws. 

Finally, we need to consider the relation of our shell 



models to the physical problems and the PDE's that mo- 
tivate these models. It is important to test the conjec- 
ture stated here in that context. In light of the above 
discussion we expect that much of what has been found 
here will translate literally to the continuous problems. 
After all, the crucial aspects are the linearity of the ad- 
vection equation, and the existence of conservation laws. 
These are unchanged in the continuous problems. Of 
course, one can expect many more numerical difficulties, 
especially due to the role of angles in the multi-point 
correlation functions. Nevertheless, the idea that the un- 
derstanding of the anomalous scaling exponent of active 
fields boils down to the analysis of eigenfunctions of a 
linear operator is expected to hold verbatim. 
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